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L  INTRODUCTION 


This  work  is  a  result  of  iirvestigations  into  methods  of  eliminating  own  ship  vibrations  from 
submarine  conformal  sonar  arrays.  Future  arr^s  may  utilize  hydrophones  that  respond  to  a  sound  wave’s 
velocity  in  place  of  traditional  pressure  sensitive  hydrophones.  In  order  to  detect  an  incoming  signal  with 
such  a  hydrophone  it  is  necessary  to  remove  signals  due  to  own  ship  vibrations  from  the  hydrophone’s 
output.  Previously  this  has  been  accomplished  using  large  metal  plates  to  damp  ship  vibrations.  Because 
of  the  high  cost  and  weight  of  these  plates,  it  is  desirable  to  instead  electronically  remove  own  ship 
vibration  signals  from  the  sensor  output. 

Because  of  changing  machinery  lineups,  ship  speed,  depth  and  other  time  varying  parameters, 
traditional  fixed  filters  are  not  adequate  to  remove  own  ship  vibration  signals  from  the  hydrophone  ouQmt. 
One  method  of  dealing  with  systems  with  time  varying  characteristics  is  with  adaptive  filter  algorithms 
which  continuously  generate  an  optimal  impulse  response  describing  the  physical  system.  Adaptive  filters 
have  been  used  in  acoustic  noise  cancellation  and  in  system  identification  for  several  years  [Ref  1,2]. 

In  this  project  both  recmsive  (UR)  and  nonrecursive  (FIR)  filters  were  used  to  predict  the  output 
of  the  following  mechanical  systems: 

1.  Torsion  waves  in  an  undamped  aluminxun  rod  (a  high  Q  system) 

2.  Torsion  waves  in  a  damped  plastic  rod  (a  low  Q  system) 

3.  Flexural  waves  in  a  damped  aluminum  rod  (a  dispersive  system) 

4.  Flexural  waves  in  a  damped  plastic  rod  (a  dispersive  system) 

5.  Longitudinal  waves  in  a  plastic  bar  (a  nondispersive  system) 

6.  Longitudinal  waves  in  a  plate  (a  high  Q,  multipath  system) 

7.  Longitudinal  waves  in  a  metal  ring  (a  high  Q,  multipath  system) 
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It  is  hoped  that  the  results  of  these  experiments  on  the  effects  of  system  physical  parameters  on 
adaptive  filter  performance  can  also  be  applied  to  more  advanced  methods  of  system  identification. 
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n.  BACKGROUND  THEORY 


A.  THE  LEAST  MEAN  SQUARE  ALGORITHM 

An  adaptive  filter  is  a  filter  that  has  the  property  of  adjusting  its  characteristics  in  response  to 
varying  ^stem  characteristics.  The  adaptive  filter  algorithm  generates  an  impulse  response,  w(t), 
describing  a  system  and  continuously  refines  it  based  on  measured  system  inputs  and  outputs.  Adaptive 
systems  have  many  applications  including  prediction,  system  identification  and  interference  cancellation. 
In  this  stutfy  adaptive  filters  utilizing  the  least  mean  square  algorithm  are  used  to  predict  the  output  of 
simple  mechanical  systems  for  a  given  input  signal.  A  block  diagram  is  shown  below; 


Mechanical  System 


System 

Input 


Figure  1  Adaptive  System  Identification  Diagram 

The  adaptive  filter  algorithm  generates  an  optimal  impulse  response  to  approximate  the  behavior 
of  a  physical  system  which  has  an  impulse  response  h(t).  In  this  stu<^  the  physical  system  is  simply  a 
vibrating  piece  of  plastic  or  aluminum  excited  by  broadband  noise.  These  systems  in  general  act  like  low 
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pass  filters  in  that  they  attenuate  high  frequency  vibrations  more  than  low  frequency  vibrations  and  are 
thus  easily  modeled  ly  adaptive  filters. 

Adaptive  filters  require  a  “desired  response”  in  order  for  the  filter  to  develop  this  optimal  impulse 
response.  In  this  case,  the  desired  response  is  the  output  of  the  physical  system,  y(t).  The  adaptive  filter 
also  utilizes  an  error  signal,  y(t)-ye(t),  where  ye(t)  is  the  adaptive  filter  output.  In  this  ystem 
identification  plication  the  least  mean  square  algorithm  is  used  to  implement  the  adaptive  filter.  The 
following  is  an  outline  of  how  this  algorithm  is  developed  [Ref  3]. 

The  error  signal  produced  in  the  aihqMive  algorithm  is 

e(t)=y(t)-ye(t) ,  (1) 

where  y(t)  is  the  ou^  of  the  physical  ystem  and  ye(t)  is  the  adaptive  filter’s  estimate  of  that  output.  The 
estimated  output,  ye(t),  is  determined  by  multiplying  the  input  x(t)  ly  the  estimated  impulse  response  w(t). 
The  vector  w(t)  is  initially  all  zeros  but  is  developed  into  an  estimate  of  the  actual  impulse  response  h(t) 
by  the  adaptive  filter  algorithm.  Using  the  tap  weight  vector  of  the  adaptive  filter,  w(t),  in  Eq.  1  results 
in; 


e(t)=  y(t)-x(t)V(t  )=y(t)-w‘(t)x(t).  (2) 

Dropping  the  time  dependence  for  clarity,  the  error  squared  is; 

e^=y^+w'xx*w-2yx‘w.  (3) 

Taking  the  e5q)ected  value  of  both  sides; 

E[e^]=E|y^]-w‘E[xxV-2E[yx‘]w.  (4) 

This  expression  can  be  simplified  by  defining  the  input  correlation  matrix,  R=E[xx‘],  and  the  cross 
correlation  matrix  P=E[yx].  Utilizing  these  definitions  in  Eq.  4  results  in  an  expression  for  the  Mean 
Square  Error  (MSE) 

MSE=E[e^]=E[y^]+w'Rw-2PSv,  (5) 
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which  shows  that  the  error  squared  is  a  quadratic  function  of  the  tap  weight  vector  w.  A  plot  of  a  MSB 
function,  known  as  an  error  performance  surface,  is  shown  in  Fig.  2  using  a  tap  weight  vector  length  of 
two.  The  MSB  is  plotted  on  the  vertical  axis  while  the  horizontal  axes  represent  possible  values  of  the  first 
and  second  elements  of  the  two  element  tap  weight  vector.  When  the  tap  weight  vector  contains  more  than 
two  elements  the  performance  surface  becomes  a  hyperparaboloid.  The  point  at  the  minimum  value  of  the 
parabolic  performance  surface  corresponds  to  the  optimal  tap  weight  vector  values.  In  this  application 
this  point  is  the  best  estimate  of  the  actual  impulse  response  of  the  physical  system  being  investigated. 


Figure  2  Error  Performance  Surface  for  Tap  Weight  Length  =  2 
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In  order  to  find  this  minimum  value  of  the  performance  surface  where  error  is  minimized,  the 
least  mean  square  (LMS)  algorithm  is  used.  Other  methods  are  available  but  the  LMS  algorithm  is  one  of 
the  most  widely  used  methods  in  adaptive  filtering  mainly  because  of  its  simplicity  and  ease  of 
computation  (the  method  does  not  require  matrix  inversion). 

The  LMS  algorithm  utilizes  the  error  defined  in  Eq.  2.  The  time  index,  t,  has  been  replaced  with 
k  to  represent  time  in  this  digital  representation  of  the  algorithm: 

e(k)=y(k)-w‘(k)x(k). 

By  using  e(k)^  fi-om  Eq.  3  as  an  estimate  of  the  mean  square  enor  in  Eq.  5  the  gradient  can  be  estimated 

by: 

Ve^=[5e/9w]^  =2e[5e/5w]=-2e*x,  (6) 


where  3e/5w=-x*  from  Eq.  2. 

This  equation  is  then  used  in  a  steepest  descent  type  of  algorithm  for  determining  the  adaptive 
filter  coefficients; 


w(k+l)=w(k)+u(-Ve^)x(k),  (7) 

where  u  is  a  positive  step  size  parameter  that  controls  the  rate  of  convergence  to  the  optimal  filter  tap 
weights.  Combining  Eq.  6  and  Eq.  7  results  in 

w(k+l)=w(k)+2ue(k)x(k).  (8) 

This  equation  is  the  basis  of  the  LMS  algorithm.  It  does  not  require  inversion,  multiplication  or 
squaring  of  a  matrix. 
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B.  INFINITE  IMPULSE  RESPONSE  AND  FINITE  IMPULSE  RESPONSE  FILTERS 


Finite  Impulse  Response  (FIR)  filters  utilize  current  and  past  values  of  the  system  input  to  predict 
system  output.  In  a  FIR  implementation  of  an  adaptive  filter,  system  output  values  are  utilized  only  as  the 
desired  output  signal.  In  Infinite  Impulse  Response  (UR)  filters,  the  tap  weight  vector  operates  on  current 
and  past  values  of  both  input  and  output  values  of  the  system  in  addition  to  using  the  system  output  as  the 
desired  signal.  This  means  that  in  an  UR  filter  in  this  application,  both  x(k)  and  y(k)  are  multiplied  the 
tap  weight  vectors  to  produce  ye(k)  where  k  is  the  discrete  time  index.  For  example  in  an  FIR  filter; 

ye(k)  =  [Xi  X2  X3  .  .  .  JtL]  [  Wi  W2  W3  . . .  Wl]‘, 

whereas  in  an  UR  filter  the  estimated  output  is  determined  from  x(k)  and  y(k); 

ye(k)=x*(k)  w(k)  +  y‘(k)  w’(fct), 

where  w’(k)  is  an  additional  weighting  vector  for  the  system  output  values.  In  a  discrete  time  system  this 
would  be  represented  by: 

ye(k)  =  [xi  X2  X3 . . .  yi  y2  ys  •  •  •  yL]‘  [  wj  W2  W3 . .  .wl  w’l  w’2  w’3  . .  .w’J. 


Stability  is  an  important  consideration  in  UR  filters  because  they  may  become  unstable  due  to  the 
presence  of  a  denominator  term  in  their  transfer  function  [Ref  8].  Thty  may  also  have  nonquadratic  error 
surfaces  or  error  siufaces  containing  local  minimums.  It  is  for  these  reasons  that  UR  filters  are  seldom 
utilized  in  conjunction  with  the  LMS  algorithm.  However,  as  will  be  discussed  in  the  results  section,  UR 
filters  offer  a  performance  improvement  without  stability  problems  for  system  identification  of  the 
nondispersive  mechanical  systems  investigated  in  this  project. 
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C.  MISADJUSTMENT 


Misadjustment  is  the  ratio  of  actual  mean  square  error  to  the  mean  square  error  that  would  be 
present  if  the  filter  coefficients  were  optimal  (at  the  bottom  of  the  parabolic  error  performance  surface). 
Because  the  LMS  algorithm  is  continually  refining  the  adaptive  filter  tap  weight  vector,  the  MSE 
fluctuates  about  the  minimum  point  on  the  error  performance  surface  after  the  filter  has  converged 
(convergence  meaning  that  the  tap  weight  vector  has  adjusted  itself  to  values  corresponding  to  the 
minimum  of  the  error  surface).  As  shown  in  reference  4,  misadjustment  is  linearly  proportional  to  both 
the  convergence  parameter  u  and  the  tap  weight  vector  length  L. 

As  u  gets  bigger,  the  LMS  algorithm  tends  to  converge  faster  due  to  the  larger  step  size.  Once 
the  algorithm  converges  on  a  minimum  error,  fiuther  iterations  cause  the  solution  to  move  erratically 
about  the  minimum  point.  A  small  u  will  minimize  the  magnitude  of  these  fluctuations  about  the 
minimum  MSE  but  will  slow  the  initial  convergence  to  this  minimum. 

As  L  increases,  the  effect  of  a  given  size  step  down  the  performance  surface  also  increases.  When 
L  is  increased  (more  terms  added  to  the  estimated  impulse  response)  the  fluctuations  around  the  minimum 
MSE  have  a  larger  effect  on  the  filter  coefficients  and  thus  increase  the  misadjustment  enor. 

D.  WAVES  IN  BARS  AND  PLATES 


This  projert  used  adaptive  filters  to  identify  the  impulse  response  of  vibrating  bars,  plates  and 
rings.  Both  dispersive  and  nondispersive  vibrations  were  studied. 

In  the  first  part  of  the  e}q5eriment,  torsion,  flexural  and  longitudinal  waves  were  selectively 
excited  electromagnetically  in  plastic  and  aluminum  rods  [Ref  5],  The  rods  were  in  a  fiee-firee  condition 
(both  ends  free  to  move).  Torsional  or  longitudinal  excitation  results  in  harmonic  modes  of  vibration 
described  by; 

fn=nc/21  (9) 
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n=mode  number=l,  2, 3... 

c=Ct  (torsional  wave  speed)  or  ci  (longitudinal  wave  speed) 
l=length  of  bar 

Flexural  waves  are  dispersive,  which  means  that  their  wave  speed  is  a  function  of  frequency. 
Free-free  boundary  conditions  were  again  used  which  resulted  in  nonharmonic  modes  of  vibration  given 
1^:  [Ref  5] 

fn=7ni^Cik/81^  (10) 

n=(3.0112,  4.9994,  7.0000,  9.0000...) 
k=radius  of  gyration  =  radius/2  for  a  round  bar 
INjar  length 

Wave  propagation  speed  is  frequency  dependent  for  dispersive  waves  with  wave  speed  given  ty: 

c=(27tfcik)'^^  (11) 

E.  THEORETICAL  OPTIMUM  IMPULSE  RESPONSE 


The  optimal  frequency  response  (which  is  the  Fourier  transform  of  the  impulse  response)  for 
single  input,  single  output  systems  is  conventionally  computed  from  [Ref  7] : 

H(f)  =  G:,y(f)/G„(f)  (12) 

Gjy  (f)  =  one  sided  cross  spectral  density 
Gxx  (f)  =  one  sided  auto  spectral  density 

As  will  be  discussed  in  the  results  section,  adaptive  filters  produced  a  better  estimate  of  the 
impulse  response  of  the  simple  systems  investigated  here  than  an  impulse  response  determined  by  taking 
the  inverse  FFT  of  the  result  of  Eq.  12. 
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in.  EXPERIMENT 


A.  ELECTRODYNAMIC  TRANSDUCTION 


Electroc^amic  transduction  was  utilized  to  generate  longitudinal,  torsional  and  flexural  waves 
in  circular  bars.  [Ref.  6] 


I . — \  PEFIMANENT  MAGNET 


COIL 


PREAMP 

AND 

HIGH 

PASS 

HLTER 


Figure  3a  Electrodtynamic  Transduction 

A  coil  of  about  20  turns  of  tliin  wire  was  wrapped  around  each  end  of  the  bar.  Each  end  of  the 
bar  was  then  placed  in  a  magnetic  field  created  by  the  two  permanent  magnets  on  each  end.  As  shown  in 
Fig.  3b  and  3  c,  changing  the  relative  orientation  between  the  permanent  magnets  and  the  coil,  different 
types  of  waves  can  be  induced  in  the  bar.  For  example,  to  induce  torsional  waves  in  the  bar  the  rod  is 
arranged  such  that  the  coil  is  lying  flat  in  the  horizontal  plane  perpendicular  to  the  magnetic  field  between 
the  two  magnets  on  either  end.  Using  the  right  hand  rule  this  will  cause  the  long  sides  of  the  coils  (which 
are  rigidly  attached  to  the  bar  with  epoxy)  to  oscillate  up  and  down  when  an  AC  current  is  passed  through 
the  coils.  [Ref  6] 
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Orientation  of  the  coil  and  magnet  pole  pieces  for  the  transduction 
01  tne  lorsiona!  mode.  The  arrows  indicate  the  magnitude  and  direction  of 
the  electromagnetic  forces  on  the  two  dominant  sections  of  the  coil  for  a 
given  phase  of  driver  current.  The  forces  produce  the  moment  that  excites 
the  torsional  oscillations. 


Orientation  of  the  coil  and  magnet  pole  pieces  for  the  transduction 
of  the  flexural  mode.  The  arrows  indicate  the  magnitude  and  direction  of 
the  electromagnetic  forces  on  the  two  dominant  sections  of  the  coil  for  a 
given  phase  of  driver  current.  The  lower  section  of  the  coil  is  in  the  stronger 
region  of  the  magnetic  field  and  hence  exerts  a  greater  force  on  the  bar  than 
the  upper  section  of  the  coil,  as  indicated  by  the  shorter  arrows  above  the 
bar.  This  leads  to  a  net  force  that  excites  the  flexural  mode. 


Figure  3b  Magnet-Coil  Orientation  for  Torsion  and  Flexural 


Waves  (From  Ref.  6) 
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Orientation  of  the  coil  and  magnet  pole  pieces  for  the  transduction 
of  the  longitudinal  (extensional)  mode.  Note  that  the  pole  pieces  have  been 
rotated  by  90°  to  concentrate  the  magnetic  field  just  beyond  the  end  of  the 
bar.  The  arrows  indicate  the  magnitude  and  direction  of  the  electromagnet¬ 
ic  forces  on  the  dominant  section  of  the  coil  for  a  given  phase  of  driver  cur¬ 
rent.  Since  this  section  of  the  coil  at  the  end  is  shorter  than  the  other  sections 
used  to  drive  the  torsional  and  flexural  modes,  it  is  useful  to  reduce  the 
amount  of  stray  magnetic  field  that  links  the  other  portions  of  the  coil  if  the 
longitudinal  mode  spectrum  is  not  to  be  “polluted”  by  the  presence  of  the 
torsional  modes. 


Figure  3  c  Magnet-Coil  Orientation  for  Longitudinal  Waves  (From  Ref  6) 


The  input  signal  into  the  input  coil  and  the  amplified  output  signal  were  recorded  using 
LabVIEW  [Ref.  11]  on  a  Macintosh  Quadra  800  computer  with  the  input  and  output  sampled  at  48  klfe. 
The  input  signal  into  the  bar  was  provided  by  the  HP  3562A  spectrum  analyzer  or  a  fimction  generator. 

Plastic  bars  (Nylon)  provided  the  best  results  for  the  purposes  of  this  proj^  because  of  its  low 
Young’s  and  Shear  moduli.  This  results  in  a  lower  sound  speed  in  the  material  which  makes  analysis 
easier  on  the  short  pieces  of  material  (less  than  a  meter  long)  that  were  used  in  this  investigation.  Copper 
and  aliuninum  bars  were  tried  but  analysis  was  difficult  because  the  waves  traveled  down  the  bars  too  fast. 

Torsion  waves  were  found  to  be  useful  for  investigating  nondispersive  waves  because  they  have  a 
lower  sound  speed  than  longitudinal  waves.  This  is  due  to  the  &ct  that  a  material’s  shear  modulus  is  in 
general  smaller  than  its  Young’s  modulus. 


B.  SHAKER  TABLE 


In  a  second  experiment,  flexural  waves  were  generated  in  a  horizontal  2.54  X  2.54  X  68  cm 
plastic  bar  displacing  one  end  of  the  bar  vertically  with  a  shaker  table.  The  input  signal  was  taken  from 
an  accelerometer  bolted  directly  on  the  shaker  table.  The  ouQmt  was  measured  with  an  accelerometer 
bolted  verticaUy  on  the  fi-ee  end  of  the  bar.  The  small  difference  in  motion  between  the  input 
accelerometer  mounted  on  the  shaker  table  and  the  end  of  the  bar  which  was  attached  to  the  shaker  table 
was  considered  to  be  part  of  the  impulse  response  to  be  modeled  the  adaptive  filter. 


C.  PIEZOELECTRIC  TRANSDUCER 


In  the  final  experiment,  two  small,  home  made  piezoelectric  transducers  were  superglued  to  the 
sides  of  a  62.4  X  15.2  X  1.5  cm  aliuninum  plate.  It  was  necessary  to  use  both  a  low  noise  preamplifier 
and  a  power  amplifier  for  the  output  because  of  the  low  sensitivity  of  these  PZT  transducers.  As  before, 
the  input  and  output  signals  were  recorded  using  LabVIEW.  One  of  the  disadvantages  of  using  the  PZT 
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transducers  is  that  th^  tend  to  attenuate  low  frequency  signals.  This  is  because  the  transducers  have  little 
of  the  mass  loading  which  would  be  required  for  them  to  operate  at  low  frequencies.  [Ref  13] 

The  accelerometers  were  also  used  to  excite  waves  in  an  aluminum  ring.  The  PZT  input  and 
ouqmt  transducers  were  attached  about  90°  apart  on  a  2  cm  thick,  10  cm  wide  aluminum  ring  with  a  17 
cm  mean  diameter.  The  ring  was  made  out  of  two  half  cylinders  which  were  bolted  together.  This 
boundary  between  the  half  cylinders  had  the  effect  of  adding  more  propagation  paths  to  the  ^stem  (the 
reflected  and  transmitted  paths  at  the  boundary).  The  ring  was  suspended  from  ring  stands  with  flexible 
cord  to  eliminate  boundary  effects  that  would  be  present  with  the  ring  laying  on  a  table. 

As  will  be  discussed  in  the  results  section,  it  was  often  necessary  to  damp  the  mechanical 
systems  to  lower  their  Q  factor.  This  was  accomplished  by  wrapping  the  material  with  vinyl  electrical  tape 
and  then  with  damping  tape  (a  heavy  metal  tape).  The  electrical  tape  provided  a  compliance  while  the 
damping  tape  provided  both  compliance  (the  sticky  part  of  the  tape)  and  mass  (the  metal  part  of  the  tape). 

The  adaptive  filters  used  to  estimate  the  impulse  responses  of  the  mechanical  systems  were 
implemented  in  Matlab  [Ref  12]  using  the  system  input  and  output  data  imported  from  LabVIEW. 
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IV.  RESULTS 


A.  ADAPTIVE  FILTER  PERFORMANCE  FOR  SYSTEMS  EXCITED  WITH  PERIODIC 

SIGNALS 


The  input  signal  used  in  most  of  the  experiments  was  white  noise.  Whenever  any  type  of  periodic 
signal  was  used  for  input  (sine  wave,  square  wave  or  sawtooth  wave)  the  adaptive  filter  predicted  the 
i^stem  output  nearly  perfectly  for  both  dispersive  and  nondispersive  systems.  Figure  4a  demonstrates  this 
for  the  case  of  a  damped  (Q=8)  plastic  bar  being  driven  with  torsion  waves  at  its  resonant  fi’equency  of  1.5 
kHz.  Figure  4b  demonstrates  excellent  performance  for  the  same  plastic  bar  being  driven  with  flexural 
square  waves  at  a  flexural  resonant  frequency  of  580  Hz. 

It  is  interesting  to  note  the  response  of  the  plastic  bar  to  the  square  wave  input.  A  square  wave 
can  be  thought  of  as  a  summation  of  sine  waves.  Because  the  bar  acts  as  a  low  pass  filter  only  the  lowest 
frequency  component  of  the  square  wave  (a  580  Hz  sine  wave)  is  present  at  the  output.  Higher 
components  are  not  likely  to  land  on  other  flexural  resonances  because  of  the  nonharmonic  nature  of 
dispersive  flexural  waves  as  shown  by  Eq.  10.  Also  note  that  in  this  example  an  adaptive  filter  tap  weight 
vector  only  24  samples  long  is  adequate  to  produce  an  estimated  ouqnit  with  low  error  even  though  the 
ouqnit  is  a  sine  wave  that  has  a  period  equivalent  to  about  82  samples. 

The  adaptive  filters  converged  (convergence  meaning  that  the  tap  weight  vector  has  a($usted 
itself  to  a  stea<fy  value  near  the  bottom  of  the  performance  surface)  in  less  than  16000  iterations  for  all  of 
the  experiments,  which  corresponds  to  0.33  seconds  at  a  48  kHz  sampling  rate. 

Convergence  times  were  determined  observing  the  adaptive  filter  behavior  for  stea^ 
sinusoidal  inputs.  This  ensured  that  the  impulse  response  of  the  physical  system  remained  the  same  during 
the  experiment.  Using  a  random  noise  input  may  have  excited  different  types  of  waves  in  the  bar  such  as 
torsion  waves  and  flexural  waves.  Figure  5a  shows  the  fi'actional  error  for  sine  wave  induced  flexural 
waves  in  a  plastic  bar  converging  in  about  3500  iterations  (equivalent  to  73  ms  at  a  48  kHz  sampling  rate) 
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Volts  .  .  .  Volts 


Input  Signal 


Actual  and  Estimated  (...)  Output  Signal 


.|q-3  Error  Between  Estimated  and  Actual  Output  Signal 


Figure  4a  1.5  kHz  Torsion  Wave  in  a  Plastic  Bar 
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Error  Between  Estimated  and  Actual  Output  Signal 


9 


to  a  steady  enor  that  is  due  primarily  to  misadjustment.  The  fractional  error  is  the  difference  between  the 
p.srimatf»d  and  actual  outputs  divided  by  the  average  rms  value  of  the  actual  output.  The  convergence  can 
also  be  seen  by  plotting  the  log  of  the  fractional  error  to  see  where  the  error  settles  on  a  steadjy  value  as  in 
Fig.  5b. 

Figure  5c  shows  the  convergence  of  a  580  Hz  torsional  square  wave  in  the  same  plastic  bar 
which  is  below  resonance  for  a  torsion  wave.  Convergence  occurs  in  about  8000  iterations  (167  ms  at  a 
48  kHz  sampling  rate).  Figure  5e  shows  the  convergence  using  the  damped  aluminum  ring  being  excited 
at  8.08  kHz  which  is  its  first  resonant  frequency.  Convergence  takes  much  longer,  at  least  250  ms, 
aiq)arently  because  of  the  multiple  propagation  paths  present  in  the  ring,  all  of  which  must  be  modeled  by 
the  tap  weight  vector.  The  high  Q  of  the  ring  (about  1000  even  after  damping)  also  slows  convergence 
because  vibrations  in  high  Q  systems  decay  slowly  resulting  in  many  passes  of  the  signal  being  present  at 
the  output  transducer.  These  two  factors  also  account  for  the  poor  performance  of  the  algorithm  in 
finding  an  impulse  refuse  that  accurately  predicts  system  output. 

The  convergence  times  for  waves  in  the  plastic  bar  can  be  compared  to  several  characteristic 
timfts  associated  with  the  bar.  The  one  way  propagation  time  for  both  torsional  and  flexural  waves  from 
the  input  coil  to  the  output  coil  is  one  of  these  characteristic  times.  Another  is  the  relaxation  time  of  the 
oscillating  bar  which  is  the  time  necessary  for  the  motion  of  the  freely  oscillating  bar  to  decay  to  1/e  of  its 
initial  amplitude.  Relaxation  time  is  computed  from  [Ref  14] 

t  =  2Q/©o  (13) 

where  ©o  is  the  natural  frequency  in  radians.  A  third  characteristic  time  is  the  “characteristic  dispersion 
time”  which  is  derived  in  the  aj^ndix.  This  is  the  time  required  for  the  second  mode  to  dephase  by  'kII 
with  respect  to  double  the  frequency  of  the  fundamental  mode  as  a  result  of  the  dispersion.  These 
characteristic  rimes  are  shown  below  for  the  0.444  m  long  plastic  bar  using  measured  values  for  torsional 
and  longitudinal  wave  speeds  of  1047  m/s  and  1425  m/s  respectively  (longitudinal  wave  speed  is  used  in 
Eq.  10  in  the  calculation  of  flexural  wave  speed). 
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Fractional  Error 


Figure  5a  Adaptive  Filter  Convergence  for  Flexural  Waves  in  a  Plastic  Bar 
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Figure  5c  Adaptive  Filter  Convergence  for  Torsion  Waves  in  a  Plastic  Bar 
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Figure  5e  Adaptive  Filter  Convergence  for  Waves  in  an  Aluminum  Ring 
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Wave  Type 

Time  to  Converge 
to  Steady  State 

Propagation  Time 

^decay  2Q/G)o 

Characteristic 
Dispersion  Time 

Torsional 

42  msec 

0.416  msec 

2.3  msec 

3.4  sec 

Flexural 

73  msec 

10.8  msec 

27.3  msec 

4.0  msec 

Table  1  Characteristic  Times  of  the  Plastic  Bar 

It  is  difficult  to  come  up  with  a  meaningful  propagation  time  for  the  aluminum  ring  because  of 
the  multiple  propagation  paths  present  in  the  ring.  However,  the  input  and  output  transducers  were  0.07  m 
apart  which  resulted  in  a  0.013  ms  propagation  time  for  longitudinal  waves  traveling  in  a  direct  path 
between  the  two  transducers. 


B.  TAP  WEIGHT  VECTOR  CONVERGENCE 


The  value  of  the  step  size  parameter  u  is  important  in  ensuring  that  the  adaptive  filter  will 
converge  to  a  stea(^  tap  weight  vector  solution  for  systems  with  time  varying  impulse  responses.  The  basis 
of  the  LMS  algorithm  is  Eq.  8: 

w(k+l)=w(k)+2ue(k)x(k). 

The  step  size  u  determines  how  big  a  jump  down  the  error  performance  surface  is  taken  during  each 
iteration.  The  jump  must  be  small  enough  to  ensure  that  the  solution  (the  tap  weight  vector  w)  can  reach 
the  bottom  of  the  parabolic  performance  surface.  (As  the  parabolic  error  performance  surface  flattens  near 
its  minimum^  big  jumps  on  the  surface  during  each  iteration  are  likely  to  miss  the  optimal  solution)  On 
the  other  hand,  u  must  be  large  enough  to  allow  the  filter  to  converge  on  a  good  tap  weight  vector 
solution  as  quickly  as  possible. 

Since  u  must  not  be  too  small  nor  too  large,  it  is  convenient  to  have  some  way  of  adjusting  u  to  a 
value  that  will  ensure  convergence  in  a  reasonable  amount  of  time  even  if  the  impulse  response  changes. 
This  was  accomplished  by  using  a  u  of  variable  size.  The  size  was  adjusted  as  a  function  of  the  error 
between  predicted  and  actual  output  with  each  iteration  [Ref  10]: 
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=  (14) 

where  p  and  y  are  constants  that  control  the  rate  of  convergence.  Note  that  u  has  units  of  inverse  power. 
This  function  is  plotted  in  Fig.  6a  using  a  variety  of  values  for  P  and  y.  A  value  of  Uo  of  0.05  was  found 
to  work  well  for  the  simple  mechanical  systems  investigated  here.  Reducing  u  prevents  the  solution  (the 
tap  weight  vector  values)  from  moving  too  for  around  the  performance  surface  minimmn  after 
convergence  and  producing  erroneous  results.  Since  this  choice  of  u  is  a  function  of  error  squared,  the  stqp 
size  will  increase  when  the  error  is  large,  such  as  when  the  physical  system  changes,  which  is  desirable  to 
aid  in  quick  adaptation  to  a  new  optimum  tap  weight  vector  solution. 

Figure  6b  demonstrates  the  ability  of  the  u  to  change  when  the  impulse  response  changes.  In  this 
case  the  transfer  function  was  altered  at  sample  2400.  At  this  point  the  ou^ut  data  record  was  multiplied 
a  factor  of  20  leaving  the  phase  unchanged  (the  actual  physical  system  was  not  altered).  Eq.  14  was 
utilized  to  adjust  u  to  quickly  find  a  new  impilse  response  to  reduce  the  squared  error.  Fig.  6c 
demonstrates  what  hrqjpened  when  a  constant  u  was  used  and  the  impulse  response  changed  at  the  2400th 
sample.  The  error  squared  grew  very  quickly  because  the  tap  weight  vector  of  the  adaptive  filter  was  not 
able  to  converge  onto  the  bottom  of  the  error  performance  surface. 

C.  TAP  WEIGHT  VECTOR  LENGTH  AND  INPUT  LAG  TIME 

The  adaptive  filter  algorithm  used  in  this  application  uses  the  system  input  x(t)  and  the  error 
signal  between  the  desired  signal  y(t)  and  the  adaptive  filter  ouqiut  ye(t)  to  determine  w(t),  the  estimated 
impulse  response  of  the  mechanical  system.  Because  propagation  of  the  vilwation  wave  firom  the  input  of 
the  mechanical  system  to  the  ouqnit  takes  a  certain  amount  of  time,  the  algorithm  is  using  information 
firom  mismatched  times.  From  Eq.  8  the  changes  in  w(t)  are  proportional  to  the  product  of  the  output 
error,  derived  from  outputs  due  to  an  input  signal  at  a  previous  time,  and  the  system  input  from  the 
present  time.  Thus  it  should  be  possible  to  delay  the  input  signal  x(t)  by: 
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Step  Size  u 
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Figure  6b  Response  to  a  Changing  Impulse  Response  for  a  Variable  u 
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Convergence  Parameter, u 

I 


Figure  6c  Response  to  a  Changing  Impulse  Response  for  a  Fixed  u 
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delay  =  Propagation  Distance  /  Wave  Speed  in  Medium 


(15) 


before  it  is  used  in  Eq.  8  to  improve  performance. 

The  propagation  time  in  each  of  the  different  systems  was  measured  using  the  HP  3 562 A 
Spectrum  Analyzer.  In  the  following  example  torsion  waves  were  electromagnetically  induced  in  a 
damped  plastic  bar  with  a  broadband  noise  input.  The  analyzer  first  determined  the  complex  frequency 
reqwnse  of  the  ^stem  and  then  took  the  inverse  FFT  of  it  to  determine  the  inqwlse  response.  The 
torsional  modes  of  the  i^stem  can  be  seen  in  the  frequency  response  on  the  far  left  of  the  plot  on 
the  bottom  of  Fig.  7.  The  iqjper  part  of  Fig.  7  shows  the  impulse  response.  Because  the  impulse  response 
is  proportional  to  the  cross  correlation  between  the  input  and  ouqiut,  peaks  in  the  impulse  response 
correspond  to  the  arrival  time  of  a  wave  at  the  output.  [Ref.  7]  The  first  pulse  rqjresents  the  propagation 
time  of  the  initial  torsion  wave  traveling  down  the  bar  and  gives  the  propagation  time  for  a  wave  traveling 
once  down  the  length  of  the  bar.  The  second  peak  represents  a  wave  traveling  down  the  bar,  reflecting  off 
the  output  end,  traveling  back  towards  the  input  end  of  the  bar  and  then  back  to  the  output  end  of  the  bar. 
The  rest  of  the  peaks  represent  further  reflections  from  the  ends  of  the  bar,  their  magnitudes  progressively 
decreasing  as  the  wave  attenuates  and  spreads  out  in  time  due  to  dispersion  (Even  though  the  waves  were 
generated  as  nominally  nondispersive  torsion  waves,  some  dispersion  occurs  because  of  the  non  ideal  end 
conditions  on  the  bar  due  to  the  imperfea  electromagnetic  coil  transducers.) 

The  first  peak  begins  at  424  psec  which  means  that  the  torsion  waves  travel  down  the  bar  in  this 
amount  of  time.  Since  the  bar  is  0.444  m  long  this  results  in  a  torsional  wave  speed  of: 

0.444  m  /  424  psec  =  1047  m/sec 

This  corresponds  to  a  theoretical  value  for  torsional  wave  speed  determined  from: 

c  =  (G/p)''" 

c  =  (0.4  X  10’°  / 1200)’^  =  1080  m/sec 
Where  G  is  the  shear  modulus  of  the  bar. 
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Figure  7  Impulse  and  Frequency  Response 


Figure  8a  shows  that  delaying  the  x(t)  input  to  the  adaptive  filter  as  suggested  Eq.  16  can 
improve  the  performance  of  the  adaptive  filter  when  modeling  torsional  wave  propagation  in  the  damped 
plastic  bar.  This  plot  shows  that  when  x(t)  is  delayed  by  about  20  samples  the  adaptive  filter  performance 
improves  for  four  different  tap  weight  vector  lengths.  The  average  error  shown  is  for  the  last  200 
iterations  of  a  4800  iteration  experiment  since  this  represents  filter  performance  after  it  has  converged  on 
a  solution.  At  the  48  kHz  sampling  rate  used  in  data  collection,  20  samples  corresponds  to  411 
microseconds,  which  is  close  to  the  propagation  time  measured  using  the  impulse  response. 

Figure  8b  shows  a  similar  result  for  a  0.92  meter  long  undamped  aluminum  rod  through  which 
torsion  waves  are  propagating.  The  aluminum  has  a  much  higher  Q  than  the  damped  plastic  bar  (600  for 
alnminiim  and  8  for  the  plastic)  which  could  be  why  the  plots  are  not  as  smooth  as  the  one  for  plastic. 
Torsional  wave  speed  in  the  bar  is  2980  m/s  which  results  in  an  optimal  delay  of  0.92m  /  2980m/s  =  0.31 
milliseconds.  This  corresponds  to  an  expected  minimum  error  at  a  lag  of  about  15  samples  at  a  48  kHz 
sampling  rate  which  is  close  to  the  results  seen  in  Fig.  8b. 

As  the  tap  weight  vector  length  L  is  increased,  the  minimum  error  decreases  as  does  the  amount 
of  lag  necessary  to  achieve  the  best  filter  performance.  This  is  because  as  the  tap  weight  length  increases 
the  input  does  not  have  to  be  delayed  as  much  for  the  tap  weight  vector  to  account  for  the  propagation 
time  of  vibrations  through  the  system.  This  is  demonstrated  by  the  “minimums”  in  the  error  squared 
stretching  out  towards  zero  lag  time  as  the  tap  weight  vertor  length  is  increased.  Each  element  of  the  tap 
weight  vector  represents  1/fs  (fs=sampling  rate)  seconds  of  time. 

Figures  9a  and  9b  show  results  for  dispersive  flexural  waves  generated  in  a  plastic  bar  for  L=128 
and  L=228  respectively.  The  minimums  occur  after  a  longer  lag  because  the  wave  speed  for  flexural 
waves  is  slower  than  for  torsional  waves  in  this  material.  The  minimum  at  a  lag  of  about  300  samples 
corresponds  to  the  travel  time  of  the  first  flexural  mode  of  the  bar  (the  slowest  of  the  flexural  modes).  This 
lag  time  can  be  predicted  by  computing  the  first  resonant  frequency  from  Eq.  10: 

f,  =  71  (2187)  (0.00318)  (3.0112)^  /  8  (0.444)^  =  126  Hz, 


34 


Avg.  Error^2  For  Last  200  Iteration: 


Avg.  Error^2  for  Last  200  Iterations 


Tap  Weight  Length  =  4  c  Tap  Weight  Length  =  8 


Figure  8b  Error  vs.  Input  Lag  for  Flexural  Waves  in  an  Undamped  Aluminum  Bar 
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Figure  9a  Flexural  Wave  in  a  Plastic  Rod  Showing  a  Minimum  at  Lag=300  Due  to  Travel 

Time  of  the  Flexural  Wave  (L=128) 
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where  ft  is  the  frequency  of  the  first  dispersive  mode  using  a  longitudinal  wave  speed  of  2187  m/s,  a 
radius  of  gyration  of  0.003 18  and  a  length  of 0.444  m.  From  Eq.  1 1  the  speed  of  the  first  mode  is: 

Cl  =  [2n  (126)  (0.003 18)(2 187)]*^  =  74  m/sec. 

This  results  in  a  travel  time  that  corresponds  to  a  lag  of  287  samples  which  is  quite  close  to  the  result  in 
Fig.  9.  The  error  with  zero  lag  is  smaller  than  the  value  that  would  be  expected  when  propagation  time  is 
taken  into  account  because  of  the  dispersive  nature  of  the  wave  which  causes  it  to  spread  out  in  time  as  it 
propagates.  This  apparently  makes  it  difficult  for  the  adaptive  filter  to  model  the  mechanical  system  with 
a  relatively  short  tap  weight  vector. 

Figure  10  shows  the  effect  of  delaying  the  input  to  the  adaptive  filter  for  the  case  of  the 
aluminum  ring,  ./^jparently,  because  of  the  multiple  propagation  paths  present  in  the  ring  and  its  high  Q 
(about  1100),  delaying  the  uq)ut  does  not  improve  filter  performance.  The  periodic  nature  of  the  result 
(about  3  samples  between  the  weak  minimums)  does,  however,  correspond  to  the  propagation  time 
between  the  input  and  output  transducers  mounted  on  the  ring.  Similar  results  were  obtained  for  the 
e?q)eriment  utilizing  PZT  transducers  moimted  on  an  aluminum  plate  excited  1^  broadband  noise. 

As  discussed  in  the  background  theory  section,  tap  weight  length  cannot  continue  to  grow 
indefinitely  to  reduce  error  because  of  misadjustment  error  which  is  inoportional  to  the  tap  weight  vector 
length.  Figure  11  demonstrates  this  for  the  case  of  broadband  noise  exciting  torsional  waves  in  the 
damped  plastic  bar.  Error  squared  initially  decreases  with  increasing  tap  weight  length  then  increases  as 
misadjustment  error  begins  to  dominate.  The  input  signal  x(t)  was  not  delayed  for  this  example.  The  first 
large  drop  in  the  error  squared  occurs  at  a  tap  weight  length  of  about  20.  The  20  tap  weights  at  a  48  kHz 
sample  rate  correspond  to  20  X  l/48kHz  =  416  microseconds  which  is  the  propagation  time  for  a 
torsional  wave  traveling  from  one  end  of  the  bar  to  the  other.  A  tap  weight  vector  more  than  20  elements 
long  is  thus  adequate  to  model  the  first  wave  reaching  the  end  of  the  bar.  The  second  large  drop  in  the 
error  squared  occurs  at  a  tap  weight  vector  length  of  about  60  which  conesponds  to  the  wave  traveling 
down  the  bar,  reflecting  off  the  end,  then  traveling  down  and  back  up  the  bar. 
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Avg.  Error'^2  for  Last  200  Samples 


Tap  Weight  Length  =  16 


Figure  10  Error  Vs,  Input  Lag  for  Waves  in  an  Undamped  Aluminum  Ring 
(Multiple  Propagation  Paths) 
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Average  Error  For  Last  200  Iterations 


Error  Vs  Tap  Weight  Length 


xio"® 
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When  periodic  inputs  were  used,  the  predicted  output  increased  drastically  as  the  tap  wei^t 
length  was  increased  past  a  certain  length  using  both  the  recursive  and  nonrecursive  filters.  It  was  evident 
that  this  effect  was  due  to  something  other  than  misadjustment  error  which  causes  enor  squared  to 
inaease  about  linearly  with  tap  weight  length  after  reaching  some  minimum  value  as  shown  in  Fig.  11. 
The  performance  of  the  adaptive  filter  improved  as  the  tap  weight  vector  length  L  got  longer  until  a 
certain  value  past  which  the  filter  was  unable  to  generate  a  stable  solution.  This  value  of  L  depends  on  u; 
u  increased  the  maximum  value  of  L  for  which  the  adaptive  filter  produced  a  stable  solution. 

An  example  of  this  phenomenon  is  shown  in  Fig.  12a.  A  580  Hz  square  wave  was  used  to  drive 
flexural  waves  in  a  damped  plastic  bar  with  the  step  size  parameter  u  initially  set  to  0.05.  It  can  be  seen 
that  inaeasing  the  tap  weight  vector  length  improved  the  performance  of  the  adaptive  filter  until  a  tap 
weight  vector  length  of  155,  above  which  the  adaptive  filter  did  not  produce  a  stable  solution.  When  u 
was  reduced  to  0.04  in  Fig.  12e-g,  the  filter  provided  a  stable  solution  when  a  tap  weight  vector  as  long  as 
166  was  used,  above  this  value  an  unstable  solution  was  again  generated. 

D.  SYSTEM  Q 

Q  is  the  quality  factor  or  “sharpness”  of  a  resonance.  It  is  equal  to  In  times  the  ratio  of  the  ener^r 
of  an  oscillator  at  resonance  to  the  energy  dissipated  per  cycle.  It  plays  an  important  role  in  the 
performance  of  the  LMS  algorithm  because  high  Q  systems  have  very  sharp  resonances  which  are 
described  by  very  long  impulse  responses.  When  these  impulse  responses  are  longer  than  the  tap  weight 
vector,  the  filter  is  insufBcient  to  model  the  physical  system  and  large  errors  between  predicted  and  actual 

ouqjut  result. 

The  Q  of  the  undamped  aluminum  bar  was  about  700  while  the  Q  of  the  undamped  aluimnum 
plate  was  about  2800  (The  Q  was  determined  by  dividing  the  resonant  frequency  by  the  width  of  the  peak 
at  its  -3  dB  down  points).  At  these  high  Q’s  the  filters  performed  quite  poorly  and  were  unable  to 
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Figure  12a  580  Hz  Sqare  Flexiu'al  Waves  in  a  Plastic  Bar  (L=24) 

Initial  u=0.05 


Volts 


4600  4620  4640  4660  4680  4700  4720  4740  4760  4780  4800 

Iteration  Number 

Figure  12b  580  Hz  Sqare  Flexural  Waves  in  a  Plastic  Bar  (L=72) 

Initial  u=0.05 


Volts 


Input  Signal 


Figure  12c  580  Hz  Sqare  Flexural  Waves  in  a  Plastic  Bar  (L=155) 

Initial  u=0.05 
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Volts 


Input  Signal 


Figure  12f  580  Hz  Sqare  Flexural  Waves  in  a  Plastic  Bar  (L=166) 

Initial  u=0.04 
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Volts 


converge  on  an  optimal  tap  weight  vector.  Filter  performance  for  the  aluminum  plate  improved  steadily 
as  the  Q  was  decreased  adding  multiple  layers  of  electrical  tape  and  damping  tape. 

High  Q  systems  such  as  the  undamped  bar  and  the  ring  (Q=l  100)  were  not  well  modeled  by  the 
adaptive  filter.  Using  longer  tap  weight  vectors  (as  big  as  2000)  and  increasing  the  number  of  iterations 
to  as  many  as  16000  produced  only  a  slight  improvement  in  performance.  However,  filter  performance 
could  be  improved  somewhat  by  reducing  the  convergence  parameter  u  by  a  factor  of  L  as  shown  in  Fig. 
13  using  data  fi'om  the  aluminum  ring.  An  optimum  value  for  L  of  248  was  used  in  this  example  by 
plotting  error  vs.  tap  weight  vector  length  as  was  done  in  Fig.  1 1 .  Reducing  u  shortened  the  jumps  down 
the  performance  surface  that  the  filter  took  during  each  iteration.  This  method  only  improved 
performance  in  the  case  of  high  Q  systems  (Q>100).  For  low  Q  systems  it  drastically  increased  the  enor  of 
the  predicted  output. 

E.  RECURSIVE  VS.  NONRECURSIVE  ADAPTIVE  FILTERS 

In  some  cases  recursive  (HR)  filters  outperformed  nonrecursive  (FIR)  filters.  Figures  14a  and  14b 
show  the  performance  of  the  adaptive  filter  (L=48)  in  predicting  the  output  for  a  low  Q  (Q=8)  system 
consisting  of  torsion  waves,  generated  firom  white  noise,  traveling  in  a  damped  plastic  bar.  It  is  clear  that 
the  HR  filter  performs  better  than  the  FIR  filter.  The  improvement  seen  when  using  the  HR  filter  is  due  to 
the  denominator  term  of  the  DR  filter  which  allows  better  modeling  of  the  resonances  present  in  the 
physical  ^stem.  Similar  results  were  obtained  for  the  other  nondispersive  systems  including  those  shown 
in  Fig.  14c  and  14d.  These  are  the  results  of  using  torsion  waves  in  a  damped  (Q=50)  aluminum  rod 
excited  broadband  noise  using  a  tap  weight  length  of  128. 

No  improvement  or  only  slight  improvement  was  seen  when  using  an  HR  filter  in  dispersive 
systems  as  seen  in  Fig.  15a  and  15b  for  flexural  waves  excited  electromagnetically  by  broadband  noise  in 
a  HampeH  plastic  bar  (L=248).  This  is  because  in  dispersive  systems  the  traveling  wave  spreads  out  in 
time  The  denominator  term  in  the  frequency  response  found  in  UR  systems  is  thus  no  longer  as  important 
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Error  Between  Estimated  and  Actual  Output  Signal 


Iteration  Number 

Figure  13a  High  Q  System,  Initial  u=0.05 
4800  Iterations,  L=248 
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Actual  and  Estimated  (...)  Output  Signal 


Volts  Volts 


Input  Signal 


Error  Between  Estimated  and  Actual  Output  Signal 
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Error  Between  Estimated  and  Actual  Output  Signal 


in  determining  the  output  from  the  input.  Figures  15c  and  15d  show  a  similar  result  for  flexural  waves  in 
a  plastic  bar  generated  tuoadband  noise  driving  a  shaker  table. 

Stability  is  an  important  consideration  in  UR  filters.  Instability  is  a  result  of  the  existence  of  poles 
in  HR  filters  which  create  a  denominator  term  in  the  frequency  response.  Performance  surfaces  of  DR 
adaptive  filters  may  have  local  minimums  and  may  even  be  nonquadratic  which  could  lead  the  tap  weight 
vector  to  converge  on  a  nonoptimal  solution.  However,  no  instability  problems  were  seen  when  using  an 
IIR  filter  on  the  simple  mechanical  systems  investigated  here.  When  instability  was  observed  (when  the 
tap  weight  vector  failed  to  converge  on  a  tap  weight  vector)  both  FIR  and  HR  filters  were  unable  to 
converge  on  a  tap  weight  vector  for  the  experiment  as  discussed  in  section  C  of  the  results. 

Both  the  FIR  and  the  IIR  filter  performed  better  on  the  data  from  the  shaker  table  experiment 
than  on  the  data  from  waves  generated  electromagnetically.  The  shaker  table  generated  “smoother” 
flexural  waves  than  did  the  waves  generated  electro^mamic  transduction.  The  coils  were  attached  to 
the  bar  1^  hand  and  could  not  be  perfectly  aligned  to  produce  purely  flexural  waves.  The  electrodynamic 
transduction  probably  produced  shear,  longitudinal  and  torsional  waves  in  addition  to  the  flexural  waves 
that  were  intended  which  made  modeling  by  the  adaptive  filter  difiicult. 

F.  DISPERSIVE  VS.  NONDISPERSIVE  PROPAGATION 


In  dispersive  systems  the  adaptive  filters  generally  behaved  poorly.  Performance  was  improved  by 
iicing  many  more  iterations  than  were  used  for  nondispersive  systems.  While  the  nondispersive  systems 
studied  required  several  hundred  iterations,  the  dispersive  systems  typically  required  at  least  16000 
iterations  to  produce  a  predicted  output  with  a  reasonable  error.  Figure  16  shows  the  result  for  flexural 
waves  induced  electromagnetically  in  a  plastic  bar  after  4800  and  16000  iterations.  Figure  17  shows  a 
similar  result  for  flexural  waves  produced  by  a  shaker  table.  IIR  filters  were  used  for  both  figures. 

An  approach  that  might  be  more  successful  in  performing  system  identification  on  dispersive 
systems  is  subband  filtering.  [Ref  9]  This  would  involve  dividing  the  input  signal  into  different  frequency 
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Figure  17b  Results  After  16000  Iterations  for  Flexural  Waves-Shaker  table  Data 

(HR  Filter) 
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bands  to  be  acted  on  an  individual  adaptive  filter  for  each  band.  A  number  of  subband  adaptive  filter 
algorithms  are  available  which  also  have  the  advantage  of  more  closely  tracking  nonstationaiy  systems. 


G.  COMPARISON  WITH  CROSS  SPECTRAL  DENSITY  METHOD  OF  SYSTEM 
IDENTIFICATION 


As  discussed  previously,  Eq.  12  is  commonly  used  to  estimate  the  impulse  re^nse  of  single 
input,  single  output  systems.  Figure.  18a  demonstrates  that  adaptive  system  identification  produces  a 
better  estimate  of  the  impulse  response  than  does  Eq.  12  for  torsion  waves  in  a  damped  plastic  bar.  The 
upper  plot  shows  the  output  produced  from  an  impulse  response  determined  an  ad^ve  filter  with  a 
128  element  tap  weight  vector.  Then,  128  element  long  samples  were  used  to  determine  Gxx  and  G^y  and 
Eq.  12  was  used  to  determine  the  firequency  re^nse.  The  inverse  FFT  of  the  frequency  response  was 
taken  to  determine  the  impulse  response,  h(t),  which  was  used  to  produce  the  middle  plot  on  Fig.  18a.  The 
bottom  plot  on  Fig.  18a  compares  the  error  associated  with  each  method.  Fig.  18b  shows  similar  results 
for  shaker  table  generated  flexural  waves  in  a  plastic  bar. 
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Actual  and  Estimated  (...)  Output  Signal-w(t) 


Error  Between  Estimated  and  Actual  Output  Signal  (...  w(t)  error,  —  h(t)  error) 


Figure  18a  Comparison  with  Cross  Spectral  Density  Method 
Damped  Plastic  Bar  in  Torsion  Excited  with  Broadband  Noise 


68 


70 


V.  CONCLUSION 


The  perfonnance  of  adaptive  filters  used  for  system  identification  of  mechanical  systems  depends 
greatly  on  the  physical  characteristics  of  that  system.  The  propagation  time,  Q  and  dispersiveness  of  the 
mechanical  system  must  be  considered  when  modeling  a  mechanical  system  with  an  adaptive  filter. 

It  is  desirable  to  take  into  account  the  time  necessary  for  vibration  waves  to  propagate  from  the 
input  of  the  system  to  the  output  so  that  the  algorithm  is  not  utilizing  data  from  mismatched  times.  This 
was  shown  to  be  quite  simple  in  the  case  of  a  torsion  wave  traveling  down  a  plastic  bar.  However,  in 
complex,  multipath  systems  determining  this  propagation  time  would  be  much  more  difficult.  In  the 
multipath  aluminum  plate  and  ring  experiments  the  performance  of  the  adaptive  filter  was  degraded  by 
delaying  the  input  to  account  for  propagation  time.  More  research  is  required  to  see  if  the  this  concqjt  is 
applicable  to  a  very  complex  system  such  as  a  submarine. 

The  high  Q  systems  investigated  were  not  well  modeled  by  the  adaptive  algorithm.  The  adai^ve 
filters  required  very  small  values  of  the  step  size  parameter  u  to  be  effective.  The  dispersive  systems 
studied  required  many  more  iterations  than  nondispersive  systems  (ly  a  factor  of  about  20)  to  determine  a 
good  estimate  of  the  mechanical  impulse  response.  Again,  more  research  is  required  to  determine  what 
role  these  factors  would  play  in  a  complex  mechanical  structure. 

The  adaptive  filter  method  of  system  identification  was  compared  to  the  more  conventional  method 
of  utilizing  the  cross  spectral  and  autospectral  densities  of  the  input  and  output  signals.  For  the  simple 
mechanical  systems  studied  here  the  adaptive  filter  performed  better. 

It  is  hoped  that  the  benefits  of  considering  the  physical  parameters  of  a  system  when  performing 
system  identification  with  adaptive  filters  will  also  prove  valuable  to  other,  more  advanced  methods  of 
system  identification. 


71 


72 


APPENDIX.  DERIVATION  OF  THE  CHARACTERISTIC  DISPERSION  TIME 


In  dispersive  ^sterns  a  pulse  will  spread  out  in  time  because  its  different  Fourier  components 
travel  at  different  speeds.  In  order  to  quantify  the  effect  of  the  dispersion  in  “stretching”  the  pulse  over 
time,  a  “characteristic  dispersion  time”  was  defined  This  constant  is  useful  in  evaluating  the  convergence 
performance  of  adaptive  filters  when  modeling  dispersive  systems.  The  derivation  is  begun  by  stating  that 
the  angular  frequency  of  the  second  flexural  mode  02  is  related  to  coi ,  the  angular  frequency  of  the  first 
flexural  mode  by: 

©2  =  2ffl,(l+82),  (Al) 

where  82  is  a  factor  to  account  for  the  anharmonic  nature  of  dispersive  waves.  Multiplying  both  sides  by 
time  results  in: 


which  can  be  reduced  to 


ffl2t  =  2©i(l+82)t, 


(A2) 


<t>2  *  2ij>i  —  2(j>)S2.  (A3) 

An  arbitrary  criteria  was  selected  to  describe  the  point  at  which  the  two  modes  were  so  out  of  phase  that 
the  pulse  was  aipreciably  spread  out  in  time.  The  time  for  this  to  occur  will  be  called  the  characteristic 
dispersion  time.  This  criteria  was  selected  to  be: 

®2t  -  2©it  =  Till.  (A4) 


Solving  for  time  gives: 


t  =  7t  /  2(©2  -  2®i). 


Since  the  angular  frequency  is  equal  to  2%  times  the  frequency  in  Hertz, 

t  =  71  /  2(27rf2  -  47rfi) 


or 


t=l/(4f2-8fi). 


(A5) 
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Computing  the  frequencies  of  the  first  two  modes  of  flexural  waves  in  the  plastic  bar  using  Eq. 
10  results  in  fi  =  81.7  Hz  and  f2  =  225.3  Hz.  From  Eq.  A5  this  results  in  a  dispersion  time  constant  of  4.0 
ms. 

The  characteristic  dispersion  time  for  the  n*  mode  to  shift  its  phase  7t/2  (a  quarter  cycle)  with 
respect  to  the  nondispersive  case  (n  times  the  fundamental  firequency)  can  be  computed  using; 

t=l/(4f„-4nf,).  (A6) 

A  characteristic  dispersion  time  was  also  calculated  for  the  case  of  torsion  wave  propagation. 
Theoretically  there  should  not  be  any  dispersion  for  torsion  waves  but  due  to  imperfections  in  the 
transducers  and  in  the  bar  itself  some  dispersion  will  occur  in  aU  practical  systems.  Additionally,  the 
derivation  of  Eq.  9  assumes  that  the  edge  of  the  bar  is  rigid  and  cannot  move  which  is  not  actually  the 
case.  The  first  two  torsional  resonant  frequencies  of  the  damped  plastic  bar  were  determined 
experimentally  to  be  1 . 1 12  kHz  and  2. 15  kHz.  Equation  A6  was  used  to  give  a  value  of  3.4  seconds  for  the 
characteristic  dispersion  time.  The  time  is  much  longer  than  for  the  flexural  wave  case  because  the  torsion 
waves  are  much  less  dispersive. 
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